Do methanethiol adsorbates on the Au(lll) surface dissociate? 
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The interaction of methanethiol molecules CH3SH with the Au(lll) surface is investigated, and 
it is found for the first time that the S-H bond remains intact when the methanethiol molecules are 
adsorbed on the regular Au(lll) surface. However, it breaks if defects are present in the Au(lll) 
surface. At low coverage, the fee region is favored for S atom adsorption, but at saturated coverage 
the adsorption energies at various sites are almost iso-energetic. The presented calculations show 
that a methanethiol layer on the regular Au(lll) surface does not dimerize. 



Self-assembled monolayers (SAM) are thin organic 
films that form spontaneously on solid surfaces. These 
systems have been the subject of intense research in re- 
cent years, both experimental and computational, see 
Ref. [J for a review. This level of interest may be as- 
cribed to their importance in wetting phenomena, tri- 
bology, chemical and biological sensing, optics and nan- 
otechnology. Among the many varieties of SAM, the ad- 
sorption process of alkanethiol molecules on the Au(lll) 
surface has been given special attention because of the 
relative simplicity of the molecules, the highly stable and 
ordered SAM structures, and the ease of preparing the 
Au(lll) surface. Despite the apparent simplicity of this 
system, its observation in various experiments has led to 
controversial results. One of the key issues is whether 
S-H bond dissociation might occur for the alkanethiol 
molecule when adsorbed on Au(lll) surface. 

Numerous experiments have focused on the alkanethiol 
adsorption on the Au(lll) surface. These efforts have 
given rise to the long-standing controversy whether 
the S-H bond of methanethiol molecules adsorbed to 
the Au(lll) surface is dissociated or not 0. Nuzzo 
et. al first observed non-dissociative adsorption of thi- 
ols on the surface by use of electron-energy-loss spec- 
troscopy (EELS), X-ray photoemission and temperature- 
programmed desorption (TPD) 2$. Based on the shape 
of the potential-energy diagram scaled with the heat 
of alkanethiol adsorption on Au(lll), it was suggested 
that S-H bond dissociation may occur for alkanethiols 
0- Kodama et al. reported thiolate radical desorp- 
tion for alkanethiol from the Au(lll) surface 4]. The 
TPD and X-ray photoelectron spectroscopy (XPS) stud- 
ies by Liu and co-workers indicated the dissociative ad- 
sorption of the methanethiol on Au(lll) 5]. The pio- 
neering theoretical work of Groenbeck and his coworkers 
showed that the S-H bond should be cleaved once the 
methanethiol molecules are adsorbed on the Au(lll) sur- 
face. In the wake of this work, it was generally assumed 
that dissociation takes place, and the cleavage of the S-H 
bond was used as a premise for further modeling of the 
methanethiol adsorbate on Au(lll) pj. However, in a re- 
cent TPD, Auger electron spectroscopy (AES), and low- 
temperature scanning tunneling microscopy (LT-STM) 
study of this system, no scission of the S-H bond Q was 



found. So far, no consistent and unified theoretical inter- 
pretation of the dissociation problem of methanethiol on 
the Au(lll) surface has been presented. Since different 
experiments, investigating the same system from differ- 
ent angles, arrived at different conclusions, theory is chal- 
lenged to propose a consistent model for the methanethiol 
adsorption process on the Au(lll) substrate. 

Guided by this motivation, in the present letter we 
propose a consistent and unified model to clarify the con- 
troversial S-H bond breaking issue. First we consider 
the adsorption energies for the non-dissociative and dis- 
sociative configurations in the c(4 x 2) superlattice at 
low coverage (0.25ML), and report for the first time that 
the non-dissociative structure is favored for the regular 
Au(lll) surface. Then we study a set of eleven non- 
dissociative structures at high coverage (1.0ML). Several 
initial configurations are designed as sulfur dimers, but 
after optimization the dimers disappear, thus we find 
that for methanethiol molecules the sulfur dimerization 
is not favored even at high coverage. To identify the 
conditions under which dissociation can occur, we make 
admission for defects in Au(lll). The S-H bond cleavage 
is found to proceed in the presence of defects in the gold 
surface. We thus report a new mechanism for the scission 
of methanethiol on Au(lll). 

Employing the VASP code y| , we follow two avenues of 
computation, using 1) the PAW potential which 
involves the generalized gradient approximation (GGA) 
[l"ll | to exchange and correlation (XC) in the framework 
of the PW91 approach; 2) the ultrasoft pseudopotential 
[I with XC based on the PBE formalism 12]. The wave 
functions are expanded in a plane wave basis with an 
energy cutoff of 400 eV. The Brillouin zone integrations 
are performed by use of the Monkhorst-Pack scheme|13|. 
We utilized a 3 x 3 x 1 k point mesh for the geometry 
optimization. The Au(lll) surface supercell consists of 
a c(4 x 2) superlattice (see Fig. QJi), i.e. 12 Au atoms 
per layer, a total of 48 Au atoms. The Au atoms in the 
top three atomic layers are allowed to relax, while the Au 
atoms in the bottom layer are fixed to simulate bulk-like 
termination|l4|-|l5j|. The vacuum region comprises seven 
atomic layers. We increased the energy cutoff to 500 eV 
and the number of k points to 8 x 8 x 1. These two cases 
differ by less than 2%. We calculated the gold lattice 
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FIG. 1: (a) The c(4x2) supercell used in the calculation, (b) 
The stable dissociative structure with the vacancy. 



TABLE I: The non-dissociative geometries and adsorption 
energies for the structures of methanethiol on Au(lll) at 
0.25ML. The entries 9, tilt direct and ds-Au refer to the polar 
angle between the normal vector of the surface and the S-C 
bond direction, the Au(lll) surface region towards which the 
S-C bond is tilted, and the shortest Au-S bond length. The 
entries initial and optimized site stand for the S atom at- 
tachment site before and after optimization. The columns 2-4 
and 6-8 list structural data pertaining to the initial and the 
final optimized geometry. The maximum adsorption energy 
is underlined. 
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constant, and found it to agree with the experimental 
value 16] within 2.1%. 

First we comment on our results related to the non- 
dissociative geometries and adsorption energies for the 
optimized configurations of the methanethiol molecule 
on the Au(lll) surface at 0.25ML, as displayed in Table 
n The units for the bond length and adsorption energy 
are Angstrom (A) and eV, respectively. Table |U shows 
that the adsorption energy for the stable non-dissociative 
structure is 0.66 eV, and the adsorption site preferred 
by the sulfur atom is located in the fee region (at fee 



TABLE II: Dissociative geometries and adsorption energies 
for the configurations at 0.25 ML. 
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but toward the atop site, top- fee), in keeping with Refs. 

Kondoh et al. observed that methanethiolate oc- 
cupies the atop site on the perfect Au(lll) surface |17| . 
but here our adsorbate is methanethiol. From Table I, it 
leans toward the atop site. The typical S-Au bond length 
is around 2.7 A, which indicates that non-dissociative ad- 
sorption may be understood as chemisorption [Isj.[7j. 

To compare dissociative and non-dissociative adsorp- 
tion energies, we include dissociative geometries and 
adsorption energies of methanethiolate on Au(lll) at 
0.25ML. Here the hydrogen atom is detached from the 
sulfur atom and forms a bond with the Au atom. As un- 
derlined in Table ITT1 the maximum adsorption energy is 
0.06 eV. The favored region for methanethiolate is the fee 
region which matches the results reported in @,0|"I13- 
In Table ITT1 there are three configurations marked by 
The first # structure for methanethiolate itself is identi- 
cal with the first one in the fee rubric, except for a differ- 
ent location of the dissociated hydrogen atom. The initial 
structures of the second and third # configuration are 
dissociative, but in the course of the geometry optimiza- 
tion, the hydrogen atom reunites with the sulfur atom. 
Thus the optimized structures are non-dissociative. 

Tableland Table ITT1 show that the adsorption energy 
of the non-dissociative structure is higher than that of 
the dissociative one by 0.6 eV. For further examination of 
our results we apply the ultrasoft pseudopotential in con- 
junction with the PBE XC to re-calculate the above 
structures, and find an adsorption energy difference of 
0.58 eV between the cases of cleavage and non-cleavage. 
Since for the perfect Au(lll) surface the non-dissociative 
structure displays a higher adsorption energy than the 
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TABLE III: The non-dissociative geometries and adsorption 
energies for the configurations at 1.0ML. 
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dissociative one by 0.6 eV, we conclude that at low or 
room temperature the adsorption of the methanethiol is 
non-dissociative, which supports the most recent experi- 
mental finding |7|. 

Next we turn to high coverage (1.0ML). We arrange 
four methanethiol molecules in the c(4 x 2) supercell 
in keeping with the experimentally detected structure 
[23 , [13 ■ The molecules labeled 1(3) and 2(4) are sym- 
metry equivalent, see Fig. The calculated results are 
shown in Table 1TTT1 The notation fcc+bri in Table llTTI dc- 
notes that the sulfur atom of methanethiol 1 is placed at 
a fee center and that of methanethiol 3 at bridge site, the 
notation hcp+fcc, top+fcc etc. is to be understood anal- 
ogously. The configurations marked by # are character- 
ized by non-dissociative initial structures but dissociative 
equilibrium structures, as emerging from geometry opti- 
mization. Their adsorption energies (-0.32 eV and 0.07 
eV) are distinctly lower than those of the alternative in- 
tact structures (around 0.20 eV), which further confirms 
our conclusion obtained for the low coverage case: the S- 
H bond is not broken when the methanethiol molecules 
interact with a perfect Au(lll) surface. At 1.0ML, the 
range of adsorption energies at various sites is within 
0.05 eV, corresponding to an almost iso-energetic situa- 
tion. These calculated results explain that the STM tip 
easily induces the motion of the methanethiol molecules 
on the surface at high coverage 0. Table ITTT1 shows that 
some distances between two S atoms are set to be less 
than 2.6 1, but after optimization, the distances become 
larger than 3.75 A. For the seventh configuration in Ta- 
ble [ffi] however, ds-s is 2.38 A. Here two of the hydro- 
gen atoms have been detached from the sulfur atoms and 
the adsorption energy is lower. Thus we conclude that 
for intact methanethiol molecules it is impossible for two 
sulfur atoms to form a dimer at high coverage, which is 
consistent with the observation reported in Refs |24j.|2fij 
and H3. 



TABLE IV: The geometries and adsorption energies for the 
non-dissociative and dissociative configurations at 0.25 ML. 
The nomenclature above or below refers to positions of the S 
or the H atom above or below the vacancy. The terms one 
bond or two bonds denote S atom bonding to one or two gold 
atoms. The term embedded means that the vacancy is filled 
by the S atom. 
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-0.14 


one bond 


2.50 


one bond 


bridge 2 Au 


2.37 
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To explore the mechanism for S-H bond breaking in 
adsorbed methanethiol, we consider an Au(lll) surface 
with defects, and for simplicity, we assume a vacancy 
in the top layer of the c(4 x 2) supercell. The calcula- 
tion shows that the formation energy of the vacancy on 
Au(lll) is 0.6 eV, and the introduction frequency of the 
vacancy is 1.5 x 10 14 Hz (2|j. The geometries and adsorp- 
tion energies for the optimized non-dissociative and dis- 
sociative configurations at 0.25 ML are indicated in Table 
IIVI It is shown that in the presence of the vacancy, the 
adsorption energy of the stable dissociative configuration 
(0.81 eV) is higher than that of the non-dissociative con- 
figuration (0.73 eV), in exact opposition to the vacancy- 
free case. The vacant site is located at the center of the 
gold atom hexagon. The most stable structure is real- 
ized when the sulfur and hydrogen atoms attach to two 
different gold atom pairs within the same hexagon, see 
Fig. If one attaches the sulfur atom on the other 

gold atoms instead of the six gold atoms surrounding the 
vacancy, the adsorption energy is almost equal to that 
found for the vacancy-free situation, since the vacancy 
exerts little effect on the S atom. The essential conclu- 
sion is that the defected Au(lll) acts as a catalyst for the 
S-H bond rupture while the perfect Au(lll) surface does 
not. It is expected that with increasing density of defects, 
the S-H bond scission becomes the prevailing mode of in- 
teraction between the methanethiol adsorbate and the 
Au(lll) substrate. Here we point out that after the loss 
of the H atom, the strength of the bond between the S 
atom and the defected Au surface exceeds that between 
the S atom and the perfect Au surface, the difference 
amounting to about 0.75 eV. This is consistent with the 
result reported in Ref. [26|. 
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FIG. 2: The PDOS projected on the H atom initially attached 
to the S atom adsorbed on defected and perfect Au(lll). 



To explore this phenomenon in greater detail, we have 
calculated the partial density of states (PDOS) for the 
H atom initially attached to the S atom, see Fig[3 For 
the perfect surface, the primary peaks of the hydrogen 
PDOS for the non-dissociative structure are at positions 
of lower energy than those for the dissociative situation. 
The electronic configuration related to the former case 
is thus more stable than that corresponding to the lat- 
ter, and the hydrogen atom forms a stronger bond with 
the sulfur atom than with the gold surface. For the de- 
fected surface, this order of stabilities is reversed. The 
primary peaks of the hydrogen PDOS for the dissocia- 
tive structure turn out to be at lower energy than those 
of the non-dissociative one, and therefore the hydrogen 
atom forms a stronger bond with the gold surface than 
with the sulfur atom. For the hydrogen PDOS related 
to the non-dissociative structure on the defected surface, 
the highest peak exceeds the second highest very sub- 
stantially. This distribution resembles more the DOS of 
the isolated hydrogen than the analogous spectrum for 



the perfect surface, indicating that the H-S bond is weak 
on the defected surface. 

In summary, we predict for the first time that the S- 
H bond remains intact when mcthancthiol is adsorbed 
on the regular Au(lll) surface, i.e., non-dissociative ad- 
sorption is thermodynamically stable. However, the bond 
breaks as admission is made for defects in the Au(lll) 
surface. At low coverage, the fee region is favored for S 
atom adsorption. At saturated coverage, the various ad- 
sorption sites turn out to be almost iso-energetic, which 
explains satisfactorily that the STM tip easily induces 
the motion of mcthancthiol molecules on the surface at 
high coverage Q. Our calculation indicates that two in- 
tact methanethiol molecules do not form a dimer on the 
perfect Au(lll) surface, which is consistent with the re- 
spective observation in Refs. [l^, and [2f|. These 
findings resolve the controversial issue of bond cleavage 
at low or room temperature. The emergence of thiolate 
on the Au(lll) surface reflects the presence of surface 
defects. Defects might be caused in the process of gold 
surface production, or temperature enhancement |24| or 
exposure to X-rays j2?j | . The sulfur dimers were proposed 
on the basis of X-ray diffraction (XRD) data obtained at 
saturated coverage [2g| , but we have found that a pair of 
intact methanethiols cannot form a dimer on the defect- 
free Au(lll) surface. We are in the process of exploring 
the conditions for dimerization, and in particular inspect 
defect and temperature related effects. 
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